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The concept of entanglement in systems where the particles are indistinguishable has been the 
subject of much recent interest and controversy. In this paper we study the notion of entanglement 
of particles introduced by Wiseman and Vaccaro [Phys. Rev. Lett. 91, 097902 (2003)] in several 
specific physical systems, including some that occur in condensed matter physics. The entanglement 
of particles is relevant when the identical particles are itinerant and so not distinguished by their 
position as in spin models. We show that entanglement of particles can behave differently to other 
approaches that have been used previously, such as entanglement of modes (occupation-number 
entanglement) and the entanglement in the two-spin reduced density matrix. We argue that the 
entanglement of particles is what could actually be measured in most experimental scenarios and 
thus its physical significance is clear. This suggests entanglement of particles may be useful in 
connecting theoretical and experimental studies of entanglement in condensed matter systems. 
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I. INTRODUCTION 

Recently there has been much interest in understand- 
ing and quantifying the entanglement present in quantum 
many-body systems. The aim of this program of research 
is to shed new light on systems, particularly strongly cor- 
related systems, that are difficult to treat with conven- 
tional approaches. Most studies have focused on quan- 
tum spin systems, especially near quantum phase tran- 
sitions, see for example 0, 0, M 01 an d references 
therein. The concept of entanglement is well-defined in 
these systems as the spins can be considered distinguish- 
able. 

The subject of entanglement becomes more subtle 
when the system to be studied consists of many indis- 
tinguishable particles. Examples of such systems span 
many fields of physics: quantum optics experiments, ul- 
tracold atomic gases, itinerant electrons and supercon- 
ductors. Even the question of which states are entan gled 
is the subject much recent debate 0, 0, E IH EH! 
E E [H [H |H|. The difficulty arises from the lack 
of individual identity of the particles that are supposed 
to be entangled, which is manifest as the necessary sym- 
metrization or antisymmetrization of the quantum wave- 
function. 

Wiseman and Vaccaro have recently proposed a 
measure of entanglement for systems of indistinguishable 
particles that is operational in the sense that it quanti- 
fies the amount of "accessible" entanglement in the sys- 
tem where a local particle number superselection rule 
restricts the possible operations that may be performed. 
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This is in contrast to other measures such as the mode 
entanglement, where the physical meaning of the entan- 
glement measure is not so clear as the measurements re- 
quired to demonstrate entanglement are not obviously 
possible. The rapidly-developing field of mesoscopic elec- 
tronics may provide a useful testing ground for compar- 
ing different notions of entanglement in condensed mat- 
ter systems as experiments to demonstrate entan glem ent 
may be feasible in the near future [Tl I20L Ell I2I I2I |24 | . 

The accessible entanglement, which is referred to as 
the "entanglement of particles" in El 1 is defined as the 
amount of entanglement that could be extracted from 
the system and placed in conventional quantum regis- 
ters, from which it could be used to perform quantum 
information processing tasks, such as teleportation. In 
many physical systems it may be difficult to extract all, 
or even some, of this entanglement, but the quantity itself 
may still give insight into the physical properties of the 
system. An analogy with thermodynamics is helpful. In 
thermodynamics it is often fruitful to consider quantities 
such as the amount of free energy in a system, even in the 
absence of an explicit scheme to extract that free energy 
If the total entanglement is taken to be analogous to the 
total internal energy, then the accessible entanglement is 
somewhat analogous to the free energy. 

In El the entanglement of particles was defined and 
evaluated for a number of states of indistinguishable par- 
ticles. However there is a lack of studies of entanglement 
of particles in explicit physical models. In this work we 
aim to fill this gap by investigating the entanglement of 
particles in a number of simple physical systems, and 
compare and contrast to other approaches to studying 
entanglement. 

We begin in Sec. ^ by reviewing the definition of 
entanglement of particles and explaining its motivation 
in terms of superselection rules and measurements. In 
Sec. IIIII we study ground and thermal states of systems 
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of bosons and fermions with a small number of modes. 
We show that for any number of modes non-interacting 
bosons have zero entanglement of particles, as one might 
intuitively expect. However non-interacting fermions can 
have non-zero entanglement of particles. We then study 
the effect of interactions on the entanglement of parti- 
cles using the Bose-Hubbard and Fermi-Hubbard mod- 
els as examples. We contrast the behaviour of the en- 
tanglement of particles with the entanglement of modes 
and show that the two measures can display opposite be- 
haviour as one varies the interaction parameter. 

In Sec. II VI we turn to multimode systems. In IIV Al we 
show how to calculate the entanglement of particles from 
correlation functions. A particularly striking example of 
entanglement of indistinguishable particles is the non- 
interacting electron gas as studied in [TH l25j . In those 
works the "two-spin reduced density matrix" — a con- 
cept common in many-body physics — is used to study 
entanglement. It is shown that there is a finite length 
over which the non-interacting electrons are entangled. 
In IIV Bl we show a similar effect in a lattice model of 
non-interacting electrons, where the entanglement per- 
sists over many lattice sites. In IIV CI we make the con- 
nection to the continuum explicit and argue that writing 
down a two-spin reduced density matrix on the lattice 
leads to difficulties in interpreting the entanglement. The 
subtle difficulty stems from the indeterminate number of 
particles at any particular location. Our results show 
that the phenomenon of entanglement of non-interacting 
electrons may feasibly be observed in an experiment. 



II. SUPERSELECTION RULES, 
MEASUREMENTS AND ACCESSIBLE 
ENTANGLEMENT 

In this section we review the concept of entanglement 
of particles, as defined by Wiseman and Vacarro [TtJ and 
explain why we consider it an appropriate measure of 
entanglement in condensed matter systems. 

In many-body physics it is common to represent the 
state of a system in the occupation-number representa- 
tion. If {ipj} is a complete set of single-particle wave- 
functions (for example modes localized in position or mo- 
mentum) then a many particle state is written 

|*> = X>|n>, w 

ft 

where ft = (ni,na, . . .) is a set of occupation numbers 
for the single particle modes (for fermions the occupa- 
tion numbers are restricted to be or 1 due to the Pauli 
exclusion principle) , and the Cf L are coefficients in the su- 
perposition. Formally, the space of occupation-number 
states is equivalent to a tensor product space where each 
mode is a factor (subsystem), and the occupation number 
of each mode represents a distinct state in that subsys- 
tem. It is thus tempting to define the "entanglement" 



in a many-body state as being with respect to this mode 
decomposition. 

Following 17], for bipartite systems we may quantify 
the entanglement of modes, Em, as 

Em(pab) = M(p AB ) (2) 

where M is some bipartite measure of entanglement (e.g. 
entanglement of formation, entanglement of distillation, 
negativity), A and B each control some subset of the total 
modes, and pab is the total state shared by A and B. 
This approach is advocated in, for example, [lti. l2(jL |27| . 
The entanglement of modes depends on which modes A 
and B control, as discussed in [2^, |2|j , but not on the 
local mode decomposition that they choose. 

The difference between entanglement of particles, 
which we define shortly, and entanglement of modes 
stems from the local particle-number superselection rule 
which may apply to systems of massive particles. Oper- 
ationally, a superselection rule (SSR) is a restriction on 
the allowed physical operations (closed or open evolution, 
preparation, measurement, etc.) on a system |30|. 

It is sometimes asserted that certain superselection 
rules apply in principle due to some underlying symme- 
try of the system, e.g. a SSR for charge that appears 
in Lorentz-invariant quantum field theories. However it 
is possible to lift superselection rules by constructing an 
appropriate reference frame for the quantity in question, 
the most famous example of this procedure being the 
thought experiment of Aharonov and Susskind [3l|. 

Whether or not superselection rules apply in principle 
is not important for our purposes. We simply note that 
often a superselection rule applies in practice due to the 
lack of an appropriate reference frame. The example we 
will be concerned with is the superselection rule for local 
particle number. Consider a bipartite state of one parti- 
cle superposed over two modes, where each party controls 
one of the modes: 

\ipe)AB = (\1)a\0)b + e^|0) A |l) s )/V2. (3) 

The phase <f) in this superposition is only meaningful rel- 
ative to some shared reference frame. Two examples 
of systems that could act as a reference frame for this 
phase are a large coherent state of light if the particle 
were a photon, or a Bose-Einstein Condensate if it were 
a bosonic atom. However without such a reference frame, 
as is generally the case in condensed matter systems, the 
phase is not accessible to experiment and the state \ij)g) 
is indistinguishable from the averaged state 

Pab = [ 7r\il>e)AB{i>e\ 
Jo 2lT 

= (\1)a(1\ ® |0) B (0| + |0) A (0| ® |1>b(1|)/2, 

which is an incoherent mixture of the particle being in 
one mode or the other. 

Notice that in the above example averaging over the 
unknowable phase <p is equivalent to projecting onto fixed 
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local particle number. This is a general result — if two 
parties, A and B, share a multimode state of indistin- 
guishable particles, Pab, and a local particle number su- 
perselection rule applies, then this state is indistinguish- 
able from the averaged state 

PAB = ^A,n B PAB^ln A ,n B , (4) 

TLA,TLB 

where TL nAi n B projects onto fixed particle number ha and 
n B at A and B, and the sum runs over all possible local 
particle numbers |T^ . 

For these reasons Wiseman and Vaccaro |T7| argue that 
the entanglement of modes does not capture the true 
amount of entanglement that the two parties, A and B, 
share since in order to take advantage of it they would 
need to be able to perform arbitrary local operations on 
the modes. In general such local operations would vio- 
late the local particle number superselection rule and are 
hence not possible in practice. 

Wiseman and Vaccaro give an operational definition 
of bipartite entanglement of indistinguishable particles 
by using the concept of a standard quantum register - 
a set of distinguishable qubits — which each party 
possesses. They define the entanglement of particles as 
the maximal amount of entanglement that the two parties 
can produce between their standard quantum registers by 
local operations on the modes that they have access to. 
Because the standard quantum registers consist of dis- 
tinguishable qubits their entanglement may be measured 
by any standard measure of bipartite entanglement. 

In |l7j only pure states are considered, however a def- 
inition of entanglement of particles Ep that applies for 
mixed states as well is 

Ep( P ab)= P nA,n B E M ( P ^' nB) ), (5) 

n A ,n B 

where Pab = ^n A .n B PAB^n A ,n B is the (unnormal- 
ized) state conditioned on obtaining the results ua and 
np for a measurement of local particle number at A and 
B, P nA ,n B = ^(Pab * ) i s * ne probability of obtain- 
ing that result and Em(p^ab ' nB ^) is the entanglement 
of modes in p( n A,n B ) gj| In words, the entanglement 
of particles is the weighted sum of the entanglement of 
modes when local particle number is measured. It is 
sensitive, for example, to entanglement in spin between 
two particles at distinct spatial locations, but not to 
"occupation-number entanglement" such as exists math- 
ematically in the state \ipg) above but would be impos- 
sible to extract in an experiment without a shared refer- 
ence frame. 

It would be most in the spirit of the operational def- 
inition to use the distillable entanglement 32] as the en- 
tanglement measure, Em- However this measure is often 
difficult to calculate, so throughout this paper we use en- 
tanglement of formation instead, as it is generally easier 
to calculate. In general the distillable entanglement is 



less than the entanglement of formation, and it is possi- 
ble for a quantum state to have non-zero entanglement 
of formation but zero distillable entanglement. 

The effective measurement of local particle number 
that appears in the definition of entanglement of particles 
is formally due to a lack of phase reference, as discussed 
above. However in many experimental scenarios the mea- 
surement does actually occur. For example in measuring 
correlations in spin between electrons in a mesoscopic 
conductor a measurement of spin (up or down) simulta- 
neously implies that an electron was also measured at 
the location of the detector. In measuring a mode or 
modes, e.g. momentum mode/s in a mesoscopic conduc- 
tor, the Hamiltonian coupling the measuring apparatus 
to the mode will typically commute with the total oc- 
cupation number for the mode/s. It is straightforward 
to show that under these circumstances the set of gen- 
eralized measurements (POVMs) [33| that can be im- 
plemented have Kraus operators that commute with the 
total occupation number. This implies that the gener- 
alized measurements obey the superselection rule. An 
example of a measurement that does not commute with 
local particle number, and therefore does not obey the 
superselection rule, is a projective measurement in the 
basis {(|0) + |l))A/2,(|0)-|l))A/2}. 

Practically it is often convenient to use entanglement 
witness to prove that a certain state is entangled. We 
note that entanglement of particles could be detected by 
measuring an entanglement witness that commutes with 
local particle number, as discussed for optical lattices 
in 34]. In other words states which are entangled in 
modes but not in particles are not detected by this type 
of entanglement witness. 



III. SOME SIMPLE SYSTEMS 

We now study some simple systems using the entan- 
glement of particles in order to build intuition for its be- 
haviour before moving to multimode systems in the next 
section. 

In order to share a state with non-zero entanglement of 
particles A and B must each be in control of at least two 
modes, and there must be at least two particles in the 
system. Therefore the simplest possible system in which 
there is entanglement of particles is two particles in four 
modes. 

In these minimal systems the entanglement of particles 
is due solely to entanglement of the modes at A and B 
when there is one particle at each location (i.e. only the 
n = 1 term from Eq. (J5J contributes). Because there are 
two modes at A and B we have an effective two qubit sys- 
tem for which it is possible to calculate the entanglement 
of formation in closed form as a function of the density 
matrix for pure |44| or mixed states |35j . 
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A. Two bosons in four modes 

Perhaps the simplest model for interacting bosons on a 
lattice is the Bose-Hubbard model. For four lattice sites 
(modes) the Hamiltonian is 



H = -t 



3 

E 



'i+i 



U^n^j-l), (6) 



where bj , bj are the usual boson annihilation and creation 

operators that satisfy [Sj,St,] = Sjji, fij = b^bj is the 
number operator for site j, and we have imposed periodic 
boundary conditions, j + 1 = for j = 3. 

For a fixed total number of bosons, N = 2 say, we write 
out the Hamiltonian matrix in the Fock basis and cal- 
culate the eigenvalues and eigenstates. Fig. ^ shows the 
entanglement of particles in the (non-degenerate) ground 
state of Eq. JBJ) as a function of U/t, where both A and 
£?'s sites are adjacent to one another. The other dis- 
tinct partition, where A and B control diagonally op- 
posite modes never contains entanglement in either the 
ground or thermal state. At U = we have no entangle- 
ment of particles, as seems reasonable since the bosons 
are non-interacting. 

In fact it is possible to show that the ground state of 
the non-interacting Bose-Hubbard model with an arbi- 
trary number of sites has zero entanglement of particles 
for any possible bi-partition. Consider an A-mode ring 
containing M non-interacting bosons. The ground state 
is 



\g(N,M)) 




VA M M! 



(7) 



where |vac) is the vacuum state containing zero particles 
in each mode. If A controls n A modes and B rip modes 
(iii + nB = A), and we project onto A and B having 
and m B bosons respectively (m A + mp = M), then using 
the commutation relations amongst the boson modes 
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oc \g{n A ,m A ))\g(n B ,m B )). (8) 

In words, the projected wavefunction is proportional to 
a factorized wavefunction where A or £Ts wavefunction 
is the ground state of their m A or rag non-interacting 
bosons as if their n A or ns modes were arranged in a 
ring. Therefore there is no entanglement of particles for 
any division of the lattice into A and B. By contrast the 
entanglement of modes is non-zero between any two par- 
titions with respect to this spatial mode decomposition, 
even for non-interacting bosons. 

For [/ ^ we have non-zero entanglement of parti- 
cles in the ground state, which increases with U jt and 



plateaus at approximately 0.1405. The limit of large U 
or small i, sometimes referred to as the hard-core bo- 
son limit, displays interesting behaviour in terms of en- 
tanglement of particles. At precisely t — the ground 
state is six-fold degenerate corresponding to the six ways 
of arranging the two bosons in four modes such that 
no mode contains two bosons. Each of these canonical 
ground states has zero entanglement of particles. How- 
ever there are linear superpositions of these ground states 
that have non-zero Ep. In particular, using degenerate 
perturbation theory we find the t — > limit of the (non- 
degenerate) ground state is the superposition 



\fj) 



1 



|1, 0,1, 0) + |0,1,0,1) --^=(|1, 1,0,0) 
+ |0,l,l,0) + |0,0,l,l) + |l,0,0,l))j/2. (9) 



The different coefficients in this sum can be understood 
as due to suppressed ability to tunnel when the particles 
are in adjacent modes; because U is much larger than 
t tunneling such that two particles end up on the same 
site is energetically unfavorable. From this expression 
we can see why the ground state has zero entanglement 
of particles in the diagonal partition (A = {1,3}, B — 
{2,4}) but non-zero Ep in the adjacent partition (A — 
{1, 2}, B — {3, 4}). The projected state for the diagonal 
partition is 

ni,i|5) oc |1,1,0,0) + |0,1,1,0) + |0,0,1,1) + |1, 0,0,1) 
= (|0,1> + |1,0)) A (|0,1) + |1,0)) B , 

which is separable. For the adjacent partition the pro- 
jected state is 

n M ls) = [|i, o,i,o) + |o,i,o,i) 

+ -±=(|0,l,l,0) + |l,0,0,l))]/2 

which is non-separable — it has Ep — h(l/2 + \/2/3), 
where h(x) — —x\og 2 (x) — (1 — x) log 2 (l — x) is the bi- 
nary entropy, and normalisation P\ A = 3/4. Therefore 
entanglement of particles is 

E P {\g)) = j/i(l/2 + V2/3) ~ 0.1405. (10) 

At non-zero temperature, T, the canonical-ensemble 
thermal state is 



p = cxp(-H/k B T)/Z, 



(11) 



where kp is Boltzman's constant. As T — > the ther- 
mal state approaches Eq. 10 for small non-zero t, and 
so should contain entanglement of particles below some 
temperature. In particular, for the ground state to have 
the majority of the weight in the thermal-state mixture 
we need kpT to be of order or less than the energy gap 
to the first excited state. 
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FIG. 1: Bipartite entanglement of particles in the ground 
state of two bosons in four modes with Bose-Hubbard Hamil- 
tonian and periodic boundary conditions. A and B control ad- 
jacent modes in the ring. Ep is zero at U — (non-interacting 
bosons) and saturates at (3/4)/i(l/2 + y/2/3) ~ 0.1405 for 
large U/t corresponding to Pi i = 3/4 and Ef = h(l/2 + 
y/2/3). 

There are four states in the t — ground-state mani- 
fold that remain at energy for small t. To first order in 
t, the state (jHJ has energy 

E g = -2V2t. (12) 

Therefore the energy gap for small t is of order 2^/2t, and 
the condition for entanglement of particles is 

k B T < 2V2t. (13) 

In Fig. |21 we plot the entanglement of particles in the 
thermal state as a function of the inverse temperature 
and the tunneling, both scaled by U to obtain dimen- 
sionless quantities. The plot is for small t/U and we see 
the type of scaling we expect — entanglement of parti- 
cles appears at inverse temperatures proportional to the 
inverse of the tunneling, as in Eq. i|13|) . 

B. Two spinless fermions in four modes 

We now study a precisely analogous model for fermions 
— non-interacting spinless fermions in a four-mode ring. 
The Hamiltonian is 

3 

# = (14) 

3=0 

where Cj, are fermion annihilation and creation op- 
erators, satisfying anticommutation relations {£j,Cj,} — 



(a) 

0.15-, 




k T/U 0.1 t/U k B T/U 01 t/U 



FIG. 2: (Color online) (a) Bipartite entanglement of particles 
in the canonical-ensemble thermal state of two bosons in four 
modes with Bose-Hubbard Hamiltonian and periodic bound- 
ary conditions. A and B each control adjacent modes in the 
ring, as in Fig. Q (b) Probability of finding one particle at A 
and one at B, P\ t \ . (c) Entanglement of formation of the a 

posteriori state, p^g ■ 

djji . Due to the Pauli exclusion principle it is not pos- 
sible to have two spinless fermions on the same site so 
there can be no on-site interaction term. 

It is straightforward to calculate the spectrum 
by Fourier transforming the annihilation operators, 
Cfc = | Ei=o e 27tijk/4 Cj , k = . . . 3. For N = 2 parti- 
cles the ground state is two-fold degenerate — a basis is 

=C\CI | vac), la,) = C^lvac)}. (15) 

As for bosons, the entanglement of particles in the ground 
or thermal state is zero for the diagonal partition. How- 
ever for the adjacent partition these basis states each 
have P M = 3/4 and E F = h(l/2 + y/2/3). Furthermore 
the equal mixture of these two states, i.e. the T — > limit 
of the canonical-ensemble thermal state, has the same 
values for Pi i and Ep 45]. So we have the somewhat 
surprising result that even non-interacting fermions can 
have non-zero entanglement of particles in the ground 
state. This is in stark contrast to bosons, where in the 
non-interacting limit the entanglement of particles was 
zero. Mathematically the reason that non-interacting 
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fermions can have non-zero entanglement of particles but 
non-interacting bosons cannot is that the commutation 
relations needed to obtain Eq. (JHJ) as a local particle num- 
ber projection from Eq. J7J do not hold for fermions. 

Motivated by this counterintuitive behaviour of non- 
interacting fermions we now turn to another simple, 
and perhaps more experimentally-relevant, model of two 
fermions in four modes — the Hubbard dimer. 



C. Hubbard Dimer 

The two-site Hubbard model (Hubbard dimer) for 
fermions with spin (e.g. electrons) is defined by the 
Hamiltonian 

H=-t J2 (4 CT Cii<r + CrJlv) + U J2 "iT^ib 
<r=T4 ' j=LM 

where j — L, R is a position label and a =f , J, is a spin 
label. The t term describes hopping between the two 
sites while conserving spin, and the U term is a coloumb 
interaction between fermions on the same site. 

The Hubbard dimer is a simple model for a num- 
ber of physical systems, including the electrons in a H 2 
molecule |33|. By varying t we have a model of bond 
breaking as the two atoms are separated. 

The ground state may be calculated exactly as a func- 
tion oiU/t, see e.g. |2q . 

| 5 ) ex Go | vac), (17) 

where 

Go = 4 T 4 +4t4i +a(£V 4 *)(4 T 4i -44 T )' ( 18 ) 



and a(x) = x + \/l + x 2 . 

For the purposes of calculating entanglement of parti- 
cles it seems most natural to imagine party A controlling 
the up and down modes of one site and party B the up 
and down modes of the other site. With this partition 
we see that the entanglement of particles in the ground 
state comes entirely from the second term in <|18fl , where 
there is one fermion at each site forming a singlet. The 
projected state, the singlet, has entanglement of forma- 
tion equal to 1, and as U/t increases the probability P\ t \ 
increases from 1/2 to 1 as the fermions are forced to lo- 
calise on each site. At a fixed temperature we see that 
the entanglement of particles reaches a peak as a func- 
tion of U/t. When interpreted as a model for H2 bond- 
breaking the peak corresponds to an optimal distance at 
which the trade-off between entanglement and probabil- 
ity of measuring one electron at each atom is maximized. 
However in reality the Hubbard model is only a good ap- 
proximation to the molecule when U ft is not too large or 
small |37j . and at finite temperature vibrational modes 
will become relevant, so it is unclear whether this effect 
could actually be observed in H^- 



The behaviour of the canonical-ensemble thermal state 
is plotted in Fig. |3| For any U/t the entanglement of 
formation of the projected state approaches 1 (singlet) 
as the system is cooled to the ground state (ksT/t — > 
0), whereas the probability, Pi,i, approaches some value 
between and 1 that increases with UJ t (as the particles 
become more localized). 

In |26j | Zanardi performed a similar calculation of en- 
tanglement in the ground state of the Hubbard dimer. 
He calculates what we refer to as the entanglement of 
modes between the two sites, which doesn't distinguish 
local entropy arising from indefinite local particle num- 
ber ("charge fluctuations") and that from entanglement 
of the spins ( "spin fluctuations" ) . From his point of view 
the ground state becomes less entangled as one increases 
U/t as it goes from a superposition over four local states 
at each site (0, f , J,, 2) to a superposition over just two 
(T,!)- From the entanglement of particles viewpoint it 
is only the "spin fluctuations" that are due to accessible 
entanglement between the two sites, and these increase 
with U/t. 

Zanardi also considers the entanglement of 
modes in the reciprocal (momentum) space, 
where the Fourier-transformed mode operators are 
Cka = {c L a + e ik *c RrT )/V2, for A; = 0,1, a =T,|. In 
this mode representation the operator that creates the 
ground state, Go, may be written 

Go = ]T [1 + e^a(U/4t)] G^G^. (19) 

k=0,l 

We can see from this expression that the two fermions are 
perfectly correlated in momentum for any U/t — both 
terms create the two fermions in the same k mode, one 
up and one down. Therefore, if A controls one k mode 
and B the other there is no entanglement of particles in 
this state. 

Finally, one could imagine A controlling both up modes 
and B both down modes. To observe this type of entan- 
glement we could, for example, use a magnetic field to 
separate up and down fermions. We could then look for 
entanglement of particles in the position or momentum 
degrees of freedom, post-selected on having one up and 
one down. From Eq. QlSj l or Eq. (|19fl we see that each 
term in the superposition has one fermion up and one 
fermion down, so Pi 1 — 1 and the entanglement of par- 
ticles coincides with the entanglement of modes. Since 
we are guaranteed to have one up and one down, the 
change between position and momentum bases is a "lo- 
cal" change of basis [46j and the entanglement is therefore 
independent of this choice. It is equal to the mode en- 
tanglement between the momentum modes, as calculated 
in |2i| 01 — ^ increases from Ep = at f//t = to 
Ep — > 1 as U/t — > 00. 

The fact that one sees entanglement in position when 
the particles are distinguished by spin or entanglement 
in spin when the particles are distinguished by position 
may be viewed as a type of "dualism of entanglement" , 
as addressed in |38| . 
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FIG. 3: (Color online) Bipartite entanglement of particles in 
the canonical-ensemble thermal state of the Hubbard dimer 
as a function of the scaled temperature, ksT/t and on-site 
interaction, U /t. A controls the up and down modes of one 
site and B the up and down modes of the other. The color 
represents the probability of finding one particle at each site 
and the height represents the entanglement of formation of 
the a posteriori state. 



Thus we see that there is rather subtle structure to the 
entanglement in the ground state of the Hubbard dimer 
that is not revealed by simply calculating the entangle- 
ment of modes. The subtlety is above and beyond the 
dependence of the mode entanglement on the choice of 
modes — it arises from considering how one might per- 
form a measurement in practice to reveal this entangle- 
ment and is captured by the entanglement of particles. 



IV. MULTIMODE SYSTEMS 

Having studied the behaviour of the entanglement 
of particles in a few small systems we now move to 
more-realistic systems containing many particles in many 
modes. This is typically the situation studied in con- 
densed matter physics; a Hamiltonian is specified in 
terms of annihilation and creation operators for bosons 
or fermions on a discrete set of lattice sites labeled by an 
index, j say. In order to begin to get a feel for the role of 
entanglement in such systems one may ask simple ques- 
tions such as: is there entanglement between the spins of 



fermions on two distinct lattice sites, ja and j'b? If we 
imagine A to have control of site ja and B to have control 
of site js then this is precisely the situation addressed by 
the entanglement of particles — A and B share a state 
of indistinguishable particles of indefinite local particle 
number. 



Entanglement of particles from correlation 
functions 



We restrict our attention to fermions for the remain- 
der of the paper as we have seen in Sec. IIIII that they 
can display counterintuitive features of entanglement of 
particles that are not seen for bosons. As an aid in an- 
swering questions such as the one posed above we show 
how to write the projected density matrices that appear 
in the definition of entanglement of particles in terms of 
correlation functions. 

First note there are four possible local states at each 
site: |0), ||), \[) and |2) corresponding to zero fermions 
on the site, a single fermion with spin up, a single fcrmion 
with spin down and a doubly occupied site. Hence if one 
traces out the rest of the lattice besides two sites one 
obtains a 16 x 16 reduced density matrix between those 
two sites. We refer to this matrix as the full two-site 
matrix in the following. 

When a local particle number superselection rule ap- 
plies, the only accessible entanglement is between the 
projected state with one and only one particle at each 
site — if even one of the sites contains either zero or two 
fermions then there is no room for entanglement as there 
is then only one possible local state at that site. This 
projected state, p^ 1 ' 1 ) is a 4 x 4 matrix that we refer to 
as the projected two-site matrix. It may calculated as 



P 



(l,i) = Tr7 



r[n M pn 



1,1 



(20) 



where p is the total state of the system (e.g. in the 
next subsection we will take p to be the grand canon- 
ical ensemble thermal density matrix), Tr ■ - g indicates 
the trace over all sites in the lattice besides ja and js, 
and n^i is the projector onto the subspace where there 
is one and only one fermion at both ja and js- The 
normalisation (trace) of p*- 1,1 -*, corresponds to the a pri- 
ori probability of detecting one and only one particle at 
each site, Pi,!. 

The projector, n^i may be written in terms of number 
operators as 

II U = n JA ^{l ~ h JAl )h jB ^(l ~ n jBl ) 
+«jat( 1 - fi jAi )n jBi {l - n jBr ) 
+n 3A i(l - n jAr )n jBr (l - n ]Bl ) 
+%4i(! - n jA] )n ]Bi {\ - n ]B] ). (21) 

The matrix elements of p^ 1 ' 1 ^ may be written as aver- 
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site element 



(TT 
(TT 
(TT 
(TT 
(U 
(U 
(U 
(XT 
(IT 
(U 



(T2 

(|2 
(12 
(2| 
(2T 
(2 1 
(22 



1 2 

|2 

2 t 
21 
2i 
22) 



correlation function 



<«jat( 1 - n 3Al) n JB^{ 1 - "jsi)) 
(n J/lT (l-n JAi )c] Bi c JBT ) 

(cLiCjAT^BTC 1 -^!)) 
("jaT( 1 - ft JAl)"-JBl( 1 ~ AjBT)) 

( c jaT C ^ iC :ibI C ^t} 

{"•ja|( 1 - "jaT^jbtC 1 - ftjsl)) 
(n J/li (l-n JAT )c] Bi c jBT ) 

(^aIC 1 - ft JAT)"-JBl( 1 - A 3bT)) 



<"-JaT( 1 - "jAl)«JBt ft JBi> 

^L| £ ^T^jbT%b|) 



{"•jaT%aI^bt( 1 



< ft JAT A iAJ.c}BT^Bl) 



(^/aT^a^jbiU 



1 jbTJ 



JA T J A i JB T jfll 



TABLE I: Above double line: the elements of the projected 
two-site matrix, p' 1 , written as correlation functions. Be- 
low double line: elements of the full two-site matrix that con- 
tribute to the two-spin reduced density matrix, but are not 
in the projected state. The elements below the diagonal are 
obtained by complex conjugation. 



ages over creation and annihilation operators as 



Tr 



3A,3b[H C j A t C j B t' C lBS"^JAS\l 



(22) 



where s, s', t, t' take the values |, I- In Table |I] these cor- 
relation functions are simplified where possible by sub- 
stituting J2U into J23)- 



B. Non-interacting electrons on a lattice 

Perhaps the simplest multimodc fermionic system one 
can imagine is non-interacting electrons in thermal equi- 
librium at zero or finite temperature. We will be inter- 
ested in the thermodynamic limit — large lattice size — 
which is of relevance to condensed matter physics, and 
serves as the starting point for more realistic, interact- 
ing, models of real materials such as superconductors. 
We aim to clarify issues of entanglement in this simple 
case in order that the same issues can be addressed in 
interacting systems. 

The Hamiltonian for non-interacting electrons is 



(23) 



<j,k> 



where < j,k > indicates that the sum runs over nearest 
neighbors j and k as defined by a link in the lattice. If we 



were to add an on-site interaction between electrons of 
opposite spin (i.e. a coloumb interaction), U J2j ^j'T^ifi 
we would have the well-studied Hubbard model, the two- 
site version of which was studied in Sec. lIII CI This sim- 
ple type of interaction might be a good starting point 
for studying how interactions affect the entanglement of 
particles. 

One's immediate reaction may be that there can be 
no entanglement in this system as the fermions are non- 
interacting and the up and down spins are indepen- 
dent. However this intuition was shown to be incorrect 
in Sec. IIII CI where we saw that the two-site version of 
this model (U = 0) contained entan glem ent of particles 
in the ground state. Furthermore, in |15ll25| it is argued 
that it is indeed possible to have "entanglement of spins 
in a non-interacting electron gas" , which is roughly the 
continuum limit of our discrete model. There are sub- 
tle differences between the entanglement as studied in 
those works and the concept of entanglement of parti- 
cles that we have focused on here. Sec. II V 01 contains a 
detailed comparison of this previous work to the current 
entanglement-of-particles approach. 

We now explicitly calculate the entanglement of parti- 
cles for non-interacting fermions on a lattice, Eq. (|23|l , in 
thermal and chemical equilibrium. We choose the simple 
case of a 1-D lattice of M sites with closed boundary con- 
ditions (i.e. a ring) for the purpose of illustration. The 
state of the system at temperature T and chemical po- 
tential p is given by the grand canonical ensemble density 
matrix 



PT = exp(-(H - pN)/k B T)/Z, 



(24) 



where N is the total number operator and 
Z is the grand canonical partition function, 
Z = Tr[exp(-(ff - pN)/k B T)\. 

In order to explicitly calculate p^ 1,1 ) we use the 
fact that the grand canonical ensemble density matrix, 
Eq. I|24(l , is a Gaussian state when the system is described 
by the non-interacting fermion Hamiltonian, Eq. 
For this reason higher-order correlation functions, as in 
Table [I] factorize into second-order correlation functions. 

Another simplifying feature is that many of the matrix 
elements are zero due to the collective SU(2) rotational 
symmetry of the model. In fact the only non-zero ele- 
ments are those along the diagonal, and the off-diagonal 

elements p^'jj = Pj/f'jx*' ^ wa y ^° see ^his is via the 
well-known result due to Weyl that states invariant un- 
der collective SU(2) rotations of the spin have the form 



p = / dUU ® UpU* ® C/t = PA U A + Ps U s (25) 



where Ha/s are the projectors onto the antisymmetric 
(spanned by the singlet) and symmetric (spanned by the 
three triplet states) subspaces, and pa/s are the weights 
of these projectors. In our case pa + Ps — Pi.i- States 
of this form are known as Werner states in quantum 
information theory. 
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The non-zero matrix elements are all determined by 
two second order correlation functions 



0.02 



^3 A ,3l 



[ n 3Al) = 
( C ] a ] C 3b1, 



2 3Al) 



1 3bV/ 



1 3bD 



(26) 
(27) 



where for the first line we have also used the transla- 
tional invariance of the lattice. The average occupation 
of any individual up or down mode which we call the fill- 
ing factor in the following, is given by n, and Cj A .j B is 
an exchange correlation between the two sites. Explicitly 
the matrix elements arc 



(l.i) 

^TT.TT 
(14) 

Pn,n 
(i.i) 



(n 2 



rt) 

2\2 



P\U1 = ( n ( 1 - n ) + \ c 3a,3b\ ) 
\„ |2 



^3 A, 3l 



C 3A ,JB I ) 

(28) 



For non-interacting fermions we can actually calculate 
the two correlation functions, n and Cj A j B , explicitly as 
a function of fi and T. The Hamiltonian (|23[) is diag- 
onal when written in terms of momentum creation and 
annihilation operators 



M-l 



H = -2tJ2 cos(2nk/M)Cl a C k(T , 



(29) 



k=0 



where 



Ck,a — 



1 



M-l 

_ e 2nijk/M c 

M f-* 

3=0 



J"- 



(30) 



The occupation of the momentum-space modes is there- 
fore 



n k - (C feT C feT ) - (C kl C kl ) - e - [Ek -ri/k B T + l 



(31) 



th 



where E k — —2tcos(2irk/M) is the energy of the k 
mode. By inverse Fourier transforming back to the 
position-space modes we obtain 



M-l 



// = 71/ 

M ^ 

k=0 



M-l 



"3A,3E 



-E 

M ^ 



^1-Ki{j A -j B )k 



n k - 



(32) 
(33) 



k=0 



Fig. 01 (a) illustrates entanglement of particles between 
two sites as the system is cooled from high temperature 
down to zero temperature (the ground state) as a func- 
tion of the inverse temperature t/ksT, for a fixed chemi- 
cal potential (which determines the T — ► filling factor) . 
At high temperature (t/ksT = 0) the probability of each 
up or down mode being occupied is 0.5 and completely 
uncorrelated with any other mode, so Pi i — 0.25 and 
there is no entanglement. As we cool to lower tempera- 
tures entanglement appears between increasingly distant 



pf 0.01 




100 



FIG. 4: (a) Entanglement of particles for non-interacting 
fermions on a 1 — D lattice with 30 sites as a function of in- 
verse temperature (t/fcsT) and separation between the sites 
(\jA — 3b\)- The subfigures (b) and (c) indicate the con- 
tributions to Ep from the probability of finding one parti- 
cle at each site, Pi.i, and the entanglement of formation in 
the a posteriori state, respectively. The chemical potential 
was /i ~ —1.89 corresponding to a T — * filling factor of 
n — 0.2. In the T — > limit we therefore roughly expect that 
Pi tl ~ (0.2) 2 = 0.04 



sites in the lattice, but the probability Pi i decreases be- 
cause the mean atom number is decreasing. The entan- 
glement in p' 1 ' 1 ) decreases with the separation between 
the two sites and goes to zero at some finite separation 
between the sites that depends on the temperature. We 
call this distance the entanglement length, r e , in anticipa- 
tion of a relationship to previous work on the free Fermi 
gas to be discussed in the next section. In this case the 
entanglement length is 5 sites in the T — > limit (ground 
state). 

Fig.[S]illustrates the entanglement of particles between 
sites in the ground state as a function of the filling factor 
n. For filling factors less than one half the entanglement 
length decreases with increasing filling factor, and by the 
time half filling is reached there is only entanglement 
between neighboring sites. For filling factors greater than 
one half the entanglement length increases again due to 
the particle-hole symmetry of the model — the filling 
factor for holes is decreasing. By contrast the probability 
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FIG. 5: (Color online) Bipartite entanglement of forma- 
tion (height) and probability, Pi,i (greyscale) as a function 
of separation and filling factor in the T — > limit of the 
grand canonical-ensemble thermal state of the Fermi-Hubbard 
model with 30 sites and periodic boundary conditions. The 
chemical potential was varied from —2 to 2 to achieve different 
filling factors. 

that one particle will be found at each site reaches a 
maximum at half-filling, as indicated by greyscale in the 
figure. 

When there are only two electrons or holes in the lat- 
tice (n = 2/(2 x 30) = 1/30) they form a singlet with 
Ep = 1 independent of the separation between the sites 
— i.e. the entanglement length is infinite. This effect is 
rather more subtle than simply entanglement between 
sites as may be seen in spin models. The wavefunc- 
tion is such that the fermions are equally likely to be 
found anywhere in the lattice (apart from on top of each 
other which is slightly more likely), but wherever they 
are found they must be in a singlet. 

Note that from the entanglement of particles perspec- 
tive entanglement is only possible if A and B each control 
both up and down modes of distinct sites. If A controlled 
the up mode and B the down mode of the same site they 
could never share any entanglement as they each only 
have one mode. This is why \ja — 3b \ = is not plotted 
in these figures. 

C. The continuum limit 

We now discuss the continuum limit of the non- 
interacting fermion lattice model and contrast the en- 
tanglement of particles with another approach that has 
been used in recent work — the so-called two-spin re- 
duced density matrix 0, [25| . 

The two-spin reduced density matrix between two 
points r and r* is defined as 

PZ*tf = (tt(r)rt(^Ar)Mr)), (34) 



where \j) s (r),ip\(r) are field annihilation/creation opera- 
tors for a particle with spin s located at position r sat- 
isfying {^ s (r),ipl,(f 1 )} = 5 SS '5(r — f 1 ). We refer to this 
matrix as the spin- correlation matrix in what follows as 
we believe the name is more appropriate. 

It was shown in 0, [2^ that for a free Fermi gas in 
thermal equilibrium the spin-correlation matrix takes the 
form of a Werner state, for the same reasons as for the lat- 
tice model in the previous section. After normalisation, 
let p — pa be the weight of the singlet. The function p 
depends on the relative distance r = |r — f*| and tem- 
perature, T; p = 1 at r = and p — > as r — > oo. The 
entanglement length, r e , is uniquely determined by 

p(r e ,T) = 1/3, (35) 

and the spins are entangled for r < r e and separable for 
r > r e . At zero temperature the relevant parameter is 
the Fermi momentum kp and the entanglement length 
scales as 

r e cx \jkp. (36) 

One might expect to see similar behaviour in a lattice 
model of non-interacting fermions, and indeed we saw in 
the previous section that the concept of an entanglement 
length persists when one considers entanglement of parti- 
cles on a lattice. We will see subsequently that the lattice 
filling factor, n, plays the role of the Fermi momentum 
in a certain limit. However first we discuss the precise 
relationship between the spin correlation matrix and the 
projected two-site matrix. 

In analogy to Eq. 134|) one may be tempted to write 
down 

Pss',W — (Cj A tCj B t'CjBs'CjAs)> (37) 

as a "two-spin reduced density matrix" in the lattice. 
However from the results of the previous section this ma- 
trix does not not correspond to the density matrix of 
a two-component quantum system in the usual sense of 
quantum information theory — the central reason be- 
ing that the two subsystems are not well-defined. By 
contrast the entanglement in Eq. llL'l'l) is experimentally 
accessible. 

To see this point first note that IL'l'l) has the same 
form as l|37|) but with projectors onto the one-particle 
subspace inserted. As we have argued in Sec. [H] the 
projector is necessary in order to define the subsystems; 
without it multi-particle correlations contribute. If one 
were to imagine extracting entangled fermions from the 
lattice then the very act of extracting is implicitly a mea- 
surement of local particle number, and entanglement can 
then exist only if precisely one fermion is found per site. 

Of course the correlation functions that make up the 
matrix could, in principle, be measured without actually 
extracting fermions |4q. Moreover, it is not difficult to 
show that entanglement in the matrix (|37|) (normalized, 
and treated as if it were a density matrix) provides a 



11 



spin element 


sum of site elements 


(TT, TT) 


(TTiTT) + (T2,T2) + (2 t,2 T) + (2 2,2 2) 


(TT.T1) 


(TUi) + (2 T,2 I) 


(TUT) 


(TUT) + (T2.12) 


(TUi) 


(TT,U) 


(TUI) 


(TUI) + (T2,T2) + (2 1,2 I) + (2 2,2 2) 


(U,1T) 


(TI, IT) 


( T i- ?J- i- ) 


(TUi) + (T2,|2) 


UUT) 


UUT) + (12,12) + (2 T,2 T) + (2 2,2 2) 


( -L T j «t -l ) 


QUI) + (2 T,2 I) 




(11,11) + (12,12) + (2 1,2 1) + (2 2,2 2) 



TABLE II: Mapping the full two-site matrix onto the spin- 
correlation matrix. The left column represents the elements of 
the spin-correlation matrix, which are obtained by summing 
elements of the full two-site matrix represented in the right 
hand column. 



lower bound on the entanglement in the projected two- 
site matrix Therefore if the matrix (|37[) were recon- 
structed experimentally, as may be possible in the near 
future in mesoscopic systems j24|, then a calculation of 
non-zero entanglement in this matrix would imply that 
the accessible entanglement would be non-zero (but the 
converse is not true). Nevertheless, it would, in our opin- 
ion, still be wrong to call the spin correlation matrix i|37|) 
a density matrix, for the reasons given above. 

There are some more intuitive reasons why we should 
not expect correlations between two fermions on one 
site and one or two fermions on the other to con- 
tribute to entanglement. The spatial wavefunctions of 
two fermions on the same site are identical and therefore 
their spin wavefunction is a singlet (as in 0, |2j| ) . By 
the monogamy of entanglement neither can be entangled 
in spin with any other. In this sense a doubly occupied 
site is like an unoccupied site — one should not include 
correlations from it in the calculation of a density matrix. 
The projected two-site matrix respects particle-hole sym- 
metry in this sense and therefore fits in naturally with 
experimental considerations in mesoscopic systems |22j . 
whereas the spin-correlation matrix does not. 

In Table [H] we show precisely the relationship between 
the two matrices by writing the elements of the spin cor- 
relation matrix as sums of elements of the full two-site 
matrix. This mapping is not a "coarse graining" in the 
sense that each element of the full two-site matrix con- 
tributes to only one element in the spin-correlation ma- 
trix — certain elements, for example (22,22), map to 
many of the spin elements. Therefore the spin-correlation 
matrix is not the "density matrix" of a well defined two- 
component system in the sense that is used in quantum 
information. We conclude that, at least in lattice models, 
the entanglement of particles — i.e. the entanglement in 
the projected two-site matrix, is what should be used in- 
stead of the entanglement in the spin correlation matrix. 

Despite these problems with using the spin-correlation 
matrix to calculate entanglement for discrete lattice mod- 



els, in the continuum limit Eq. i|34[l recovers a interpreta- 
tion as a density matrix of two spins. The reason for this 
is basically that the probability of finding two electrons 
at a particular location in space is negligible compared 
with the probability of finding one electron. 

To see this in more detail, let N non-interacting elec- 
trons be confined to a region [0, L] in one dimension. De- 
fine a set of M orthonormal wavefunctions, {ipj (x) , j — 
... M — 1}, on the region by 



l/Ve, x £ [je, 
0, x £ [je, (j -\ 



(i + l)e] 

•1)6] 



(38) 



where e = L/M. Let Cj a be the annihilation operator 
that destroys an electron in the jth region with spin a. 
For non-interacting electrons in thermal equilibrium the 
probability of finding one electron in the jth region is 
{fija) = Ne/L, while the probability of finding two elec- 
trons, one up and one down, scales as {fij^hji) = 0(e 2 ). 
Therefore in the e — > limit, where ipj{x) approaches a 
delta function at x = je, the probability of finding two 
electrons at the same site becomes negligible compared 
to the probability of finding one. Hence the projector 
in Eq. 12211 has no effect, and the projected two-site ma- 
trix approaches the spin-correlation matrix (|34l) in the 
continuum limit. 

The question of using l|34|l as a density matrix was 
considered in the appendix of [39|. Their explanation 
agrees with ours for the case where there is one and only 
one particle in each of local modes. Our results show 
why this is a valid assumption in the continuum limit, 
but also apply to more general situations when there is a 
non-zero probability of finding two electrons in the same 
local mode. 

With this calculation in mind we should expect that 
in the limit of low filling factor on the lattice the differ- 
ence between Eq. (I22|l and Eq. I|37|l will be negligible. 
We define the entanglement length for the lattice as the 
smallest | ja — 3b | for which the entanglement of the two- 
site matrix is zero. Fig.[fj]shows the entanglement length 
versus the inverse of the filling factor alongside the "en- 
tanglement" length calculated from the spin correlation 
matrix. We see that in the limit of small filling factor 
(l/n large) the entanglement length scales linearly with 
1/fi, as one might expect from 0,113 since kp cx n in 
one dimension. 

For non-interacting electrons we see from Tab. ITTI that 
the many-electron correlations simply add a term propor- 
tional to the identity to the density matrix. This can only 
have the effect of diluting the entanglement and decreas- 
ing the entanglement length. We see from the inset in 
the figure that the projected matrix therefore predicts a 
longer entanglement length than for the spin-correlation 
matrix for some values of the filling factor approaching 
half-filling (n = 1/2). For 0.45 < n < 0.55 (not plot- 
ted in the figure) the spin-correlation matrix would have 
predicted an entanglement length of 1 site (i.e. not even 
nearest- neighbors are entangled), whereas the projected 
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(filling factor) 

FIG. 6: Entanglement length (crosses), r e , versus inverse fill- 
ing factor, 1/n, for non- interacting fermions on a 1 — D lat- 
tice. The circles show what the entanglement length would 
have been if it were defined using the spin-correlation matrix, 
Eq. 13711 . instead of the projected matrix, Eq. (1221 . The line is 
a linear least-squares fit to the crosses. The inset shows that 
the projected matrix predicts a longer entanglement length 
for some values of the filling factor. 

matrix gives r e = 2 (nearest- neighbors are entangled). 
This is precisely the limit where there is a good chance 
of finding two electrons at the same site. 

V. SUMMARY AND CONCLUSION 

An idea that has attracted much attention recently 
from the quantum information theory community is that 
the improved understanding of entanglement that has 
been developed may lead to new insights into the physics 
of strongly-correlated systems in condensed matter. Typ- 
ically in these systems the particles are itinerant and so 
must be treated as indistinguishable. Thus if one wishes 
to explore the role of entanglement one must have a good 
understanding of what it means for indistinguishable par- 
ticles to be entangled. 

We have argued that in many situations of interest 
in condensed matter a local particle number superselec- 
tion rule applies to operations, such as measurements, 
that can be performed on the system. Thus the notion 
of entanglement of particles, introduced in |17| . may be 
a more appropriate measure of entanglement to use in 
studying these systems than the entanglement of modes 
which appears in many previous studies. The physical 
meaning of entanglement of particles is clear — the sub- 
systems that are supposed to be entangled are estab- 
lished by measurement of local particle number by the 
two parties. The mathematical entanglement in occupa- 



tion number (mode entanglement) may still display in- 
teresting behaviour, however its physical significance is 
less clear as the measurements that could be performed 
to observe it are unspecified. 

In order to get a feel for how the entanglement of parti- 
cles compares with the entanglement of modes we began 
by studying some simple systems that are analytically 
solvable. The minimal situation in which entanglement 
of particles is possible is two particles in four modes. 
We found that the entanglement of particles was zero for 
bosons but non-zero for fermions for two non-interacting 
particles in a four-mode ring. In both cases the mode 
entanglement was non-zero. For the Hubbard dimer we 
calculated the entanglement of particles according to a 
number of different mode decompositions and compared 
with previous studies of mode entanglement |26|. 

Finally we studied non-interacting fermions on a lat- 
tice and compared with previous results for this system. 
We first showed how to write the projected matrix for 
one fermion on each of two distinct sites in terms of cor- 
relation functions. In agreement with previous results 
regarding the free Fermi gas 0, |25| we found an "entan- 
glement length" in the system beyond which the fermions 
are not entangled. It is intriguing that this length extends 
over multiple lattice sites and persists even when one con- 
siders the more restrictive criteria of entanglement of par- 
ticles where real measurements are considered. Thus the 
phenomena of entanglement of non-interacting fermions 
should be experimentally observable, perhaps in opti- 
cal lattice set-ups where condensed matter Hamiltonians 
may be engineered. It cannot be dismissed as trivially 
due to the antisymmeterization of the wavefunction and 
unobservable. Finally we showed precisely how the en- 
tanglement of particles relates to the two-spin spin re- 
duced density matrix |l5ll25j| in the continuum limit. 

Recently there has been interest in studying scal- 
ing laws for entanglement entropy in arbitrary dimen- 
sions 0, 0| . The entanglement entropy corresponds to 
the entanglement of modes for pure states when entangle- 
ment of formation is used as the measure. An interesting 
direction for future research would be to see if such scal- 
ing laws persist when the more-restrictive criteria of en- 
tanglement of particles is used, rather than entanglement 
of modes. 

In conclusion, we believe that the entanglement of par- 
ticles may be a useful concept to consider alongside mode 
entanglement in studying systems of indistinguishable 
particles that are central to condensed matter physics. 
We have showed that it is physically well-motivated by 
measurement considerations and leads to distinct phe- 
nomenology of entanglement in a few simple systems. 
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